Direct exchange in the edge-sharing spin-1/2 chain compound MgVOg 
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Bandstructure calculations with different spin arrangement for the spin-chain compound MgVOs 
have been performed, and paramagnetic as well as magnetic solutions with ferro- and antiferromag- 
netically ordered chains are found, the magnetic solutions being by 0.22 eV per formula unit lower 
than the paramagnetic one. The orbital analysis of the narrow band crossing the Fermi level in the 
paramagnetic solution reveals that the band has almost pure vanadium 3d character, the lobes of 
the relevant d-orbitals at the neighboring in-chain sites being directed towards each other, which 
suggests direct exchange. The tight-binding analysis of the band confirms the strong exchange trans- 
fer between neighboring in-chain V-ions. Besides, some additional superexchange transfer terms are 
found, which give rise both to in-plane coupling between the chains and to frustration, the dominant 
frustration occurring due to the interchain interactions. 
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Spin chains and ladders are of fundamental intere^ 
for solid state physics due to their peculiar properties.El 
In the last yearSptaany spin-chain ccmipounds were found, 
mostly cupratesErtl and vanadatesD There one can dis- 
tinguish between the corner-sharing compounds withj^ 
180° (T)ransition metal-(L)igand-T bond (SrsCuOjfl) 
and the edge-sharingxompounjds with a 90° T-L-T bond 
(LisCuOs, CuGeOall MgYoJ). In cuprates the relevant 
d orbitals are directed towards the ligand ions, which re- 
sults in a strong antiferromagnetic superexchange inter- 
action for a 180° T-L~T bond and a weaker ferromag- 
netic coupling for a 90° T-L-T bond according to the 
Goodenough-Kanamori-Anderson (GKA) rules.l3 Devia- 
tions from the GKA rules are known for CuGeOs due to 
the presence of side groups .13 

Recently, the compound MgVOs was proposed as a 
new candidate for a model spin-chain system ancLmag- 
netic susceptibility measurements were presented^ The 
data suggest short range antiferromagnetic spin correla- 
tions with the constant of the high-temperature Curie- 
Weiss law 6 sa —100 K. The data were analyzed within a 
ID spin-1/2 Heisenberg model with the nearest neighbor 
Ji and the next nearest neighbor J2 exchange coupling 
and a frustration a = J2I J\ close to the critical valueEl 
ac — 0.24 was found. Here we present bandstructure cal- 
culations for this compound and a corresponding orbital 
and tight-binding (TB) analysis. 

The base-centered orthorhombic crystal structure of 
MgVOa was determined by Bouloux et aln and is shown 
in Fig. [|. It consists of edge-sharing VO2 chains running 
along the y-direction which are coupled in x-direction 
by V-O-O-V bonds. For a convenient presentation of 
the orbital analysis the coordinate system is. rotated 90° 
about the x-axis against the standard onetij so that the 
space group reads as Bm2ib instead of Cmc2i given in 
Ref. |. 
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FIG. 1. Crystal structure of MgVOa. The dashed lines 
depict the Bravais unit cell. The lattice parameters of the 
cell are a = 5.243 A, 6 = 5.293 A, c = 10.028 A. 

One notes a slight tilting of the VO5 pyramids out 
of the z-direction and an asymmetric coordination of V 
and 01 ions. If one distorts the structure in a way de- 
picted in Table || both the tilting and the asymmetry are 
removed, a center of inversion appears, and the space 
group becomes Bmmb. All the following results were ob- 
tained using the model structure, because this simplifica- 
tion considerably reduces the computational efforts, and, 
as we have checked for the non-spin-polarized case, the 
deviation from the result obtained with the real crystal 
structure is negligible. 



TABLE I. The Wyckoff positions of atoms in MgVOs 
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The calculation of the bandstructure was per- 
formed using the full-potential nonorthogonal local- 
orbital minimum-basis bandstructure scheme (FPLQ^ 
within the local spin density approximation (LSDA).O 
The calculation was non-relativistic, the exchange and 
correlation potential was taken from Ref. |lj. The set 
of valence orbitals was chosen to be Mg: 3s3p3(i, V: 
3s3p4s4p3(i and O: 2s2p'id. The inclusion of the vana- 
dium 3s3p states turned out to be unavoidable since the 
V-02 distance of about 3.11 ao is small enough to yield a 
slight overlap of these states. The oxygen and magnesium 
id orbitals were taken to increase the basis completeness; 
though being not occupied they contribute to the overlap 
density. The extent of the basis orbitals, controlled by a 
confining potential (r/ro)"*, was optimized with respect 
to the total energy. 

In our calculation the Bloch state \\<iv) is composed 
of overlapping atomiclike orbitals \Lij) centered at the 
atomic sites j in the unit cell i with coordinates Ry : 
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with the normalization condition (ki/|kz^) — \. L stands 
for {nlma} denoting the main, orbital, magnetic quan- 
tum numbers and the projection of spin, respectively. 

We have performed three kinds of computation: one 
non-spin-polarized and two spin-polarized with coUinear 
and anticollinear spin polarization at the neighboring 
in-chain vanadium ions. In the last case we assumed 
ferromagnetic order along the x- and z-directions. In 
both spin-polarized calculations we have found magnetic 
(ferro- and antiferro-) solutions with the total energy be- 
ing about 0.22 eV per formula unit lower than that of the 
paramagnetic solution of the non-spin-polarized calcula- 
tion. The LSDA accuracy does not allow to determine 
which of the magnetic states is preferable. The magnetic 
moment (n^ — n|) /iB of the vanadium ion in both mag- 
netic solutions is close to the saturated value 1 /^b • The 
magnetic moments of the other ions are much smaller, the 
only appreciable one occurs at apex oxygen (02) having 
a value of about 0.05 /ib and being directed opposite to 
the moment of the neighboring vanadium ion. 

The results of all three calculations are presented in 
Fig. 0. The paramagnetic solution has metallic charac- 
ter with a half-filled conduction band at the Fermi level, 
whereas the band splits in two in the magnetic solu- 
tions and an insulator gap opens, being about 0.5 eV and 




FIG. 2. Band structure and density of states of MgVOs 
for the paramagnetic (a), ferromagnetic (b) and antiferro- 
magnetic (c) solutions. The symmetry point coordinates of 
the B-base centered orthorhombic Brillouin zone are given in 
units 7r(2/a, 1/6, 2/c). The majority and minority (tinier dots 
and shaded DOS) spin parts are shown. 



0.8 eV in ferro- and antiferromagnetic cases, respectively. 
However, it can be expected that the real gap is caused 
mainly by the electron correlation and is considerably 
larger than that one produced by the magnetic ordering. 
The correlation gap should persist in the paramagnetic 
case as well. 

For a deeper understanding of the electronic structure 
we have performed an orbital analysis for the paramag- 
netic case. The weight W^'^ of the orbital \Lj) in the 
Bloch state |kj/) (Eq. |l|) was taken as 



w^'; 



Eict^ 



Lij\ 



(2) 



The sum of all weights X^l i ^Lj i^ approximately unity, 
with some deviation due to the nonorthogonality of the 
basis orbitals. Actually the orbital weights were normal- 
ized with respect to the sum. 
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FIG. 3. Dispersion (a) and orbital weights of the upper (b) 
and lower (c) parts of the band crossing the Fermi level in the 
paramagnetic solution. The coordinates of fc-points are given 
in the same units as in Fig. ti. 

According to the analysis, the valence bands consist 
of a lower oxygen 2p-orbital complex which is separated 



by 3 eV from upper vanadium Sd-orbitals (cf. Fig. ||). 
The 3(i-orbitals are split according to the standard crys- 
tal field rules. At (0,1,0) their energy rises in the order: 
d.j.2_y2^ dyz, dzx, d^z'^-r^i dxy Only the lowest Zd^^^y-^- 
orbital is partially occupied, being half-filled with two 
electrons per two vanadium atoms in the primitive unit 
cell. The corresponding band (see Fig. H) is very nar- 
row having a width of only 0.8 eV. The weight of the 
V: 3dx2-y2 orbital in the band is larger than 70 per- 
cent, but there are considerable contributions up to 20 
percent from the 01: 2px and the V: 3dyz orbitals, as 
well as smaller contributions from some other orbitals. 
Comparing the situation with that in cuprates having 
edge-sharing Cu02 chains (Li2Cu02 or CuGeOs) one 
observes an important distinction: in the cuprates the 
relevant Cu: 3d-orbital is directed towards the oxygen 
ions, whereas in the present case it is directed towards 
the neighboring vanadium ions. It implies an indirect 
superexchange mechanism via the intermediate ligand in 
the cuprates but probably a dominant direct exchange 
process between neighboring vanadium ions and a much 
smaller superexchange hopping to the second in-chain 
neighbor in MgVOa. Let us note that the different co- 
ordination of the relevant Bd-orbital in vanadates and 
cuprates with similar crystal structure has to be a rather 
common feature, because in vanadates the energetically 
lowest Bd-orbital is half-filled in contrast to the highest 
one in cuprates. 

The coupling between the neighboring chains occurs 
mainly in the xy-plane via the intermediate 01: 2px or- 
bitals. This is manifested by the strong dispersion along 
the {kx, 0, 0) direction (see Fig. ||) of the lower part of the 
band, which has a remarkable contribution from the or- 
bital, whereas the upper part having a contribution only 
from the 01: 2py orbital is almost dispersionless. Fig. 
showing the relevant orbitals at the F-point illustrates 
the interchain coupling via cr-bonds of neighboring 01: 
2px orbitals. 

To quantify the result a TB analysis has been per- 
formed for the relevant band of the paramagnetic solu- 
tion. The hopping processes to the nearest (iij,) and the 
next nearest (^2^) in-chain neighbors, as well as two hop- 
ping processes {tx,txy) to the next chain and one (txyz) 
possible coupling in the z-direction were taken into ac- 
count: 



Eq — -Ek = 2tiy cos ky — h 2t2y cos kyb + 2tx cos kxa 

, b a b c 

itxy cos kx a cos ky — \- otxyz cos kx — cos ky — cos kz — , 



yielding the parameters (in meV): t 



ly 



125, t 



2y 



20 ... 26, ^a; = 50, txy = 20 and a tiny value of 3 for txyz- 
The data confirm that the ratio t2y/tiy, a measure for the 
frustration in the chain direction, is much smaller ths-a, 
the corresponding value for CuGeOa and Li2Cu02ii3'O 
It should be noted that the parameter tiy is much larger 
than the corresponding inter- ladder hopping in NaV205 



given in Ref. HS, probably due to a mutual compensation 
of different transfer paths in the ladder compound. 

An estimate of the corresponding exchange integrals of 
the effective Heisenberg Hamiltonian 



H — - } JjSjSi+j 



can be performed in the one-band Hubbard description, 
which gives an antiferromagnetic exchange of roughly 
Jj = 4:t'^/U with a value of U still to be determined. 
Using the experimental value of Jiy sa 10 meV given in 
Ref. ^ one gets U = 6.25 eV which seems to be slightly 
overestimated. Most probably, additional ferromagnetic 
processes (possibly via the Sdyz orbital or indirect pro- 
cesses via oxygen) have to be included to improve the 
estimate. Our TB analysis suggests a rather large value 
for Jx- It gives two-dimensional antiferromagnetic or- 
der with the frustration terms J2y and Jr^y (see Fig. 0). 
Though t2y and t^y are of nearly the same value, due 
to the larger coordination number, J^y dominates. In a 
simple mean field approach one can define an effective 
frustration exchange Jeff = 'ZJxy + J2y yielding an effec- 
tive frustration Jes/Jiy of roughly 0.1 which is a very 
preliminary analysis, however. 



transition metal (V) ions. It suggests a direct antiferro- 
magnetic exchange process between the neighboring in- 
chain spins. A similar process was proposed eatlier for 
the inter-ladder exchange coupling in MgV205£3 The 
importance of the direct vanadium d-d transfer was also 
anticipated in Ref. ^for CaV409, a compound with V-V 
distance slightly larger (3.00 A) than in MgVOa (2.96 A), 
where a value of 80 meV for the transfer was reported. 
Besides the direct exchange we have found also some ad- 
ditional superexchange terms which give rise to the cou- 
pling between the chains in the x-direction and to the 
frustration. The dominant frustration occurs between 
the neighboring chains in contrast to the naive picture of 
the frustration due to the next-nearest in-chain neighbor 
superexchange. It suggests that the experimental data 
require an analysis in the framework of a two- rather 
than a one-dimensional Heisenberg model. 



ACKNOWLED GMENTS 



We thank Christoph Geibel for useful discussions. 




FIG. 4. Relevant V: 3d^2_y2 and Ol: 2px,y orbitals cou- 
pling the spins in chain j/-direction and in the orthogonal 
a;-direction and the chosen in-plane transfer paths. The or- 
bital phases correspond to the F-point. 
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Thus, we have found that the VO2 chains of MgVOs 
have indeed a spin-1/2 structure. In contrast to the crys- 
tallographically similar edge-sharing cuprates the rele- 
vant 3d orbital is directed towards the nearest in-chain 



